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ABSTRACT 

Red optical and near-infrared spectroscopy are presented for SDSS J125637. 13—022452.4, one of 
only four L subdwarfs reported to date. These data confirm the low-temperature, metal-poor nature 
of this source, as indicated by prominent metal-hydride bands, alkali lines, and collision-induced H2 
absorption. The optical and near-infrared spectra of SDSS J1256— 0224 are similar to those of the 
sdL4 2MASS J16262034-f 3925190, and we derive a classification of sdL3.5 based on the preliminary 
scheme of Burgasser, Cruz, & Kirkpatrick. The kinematics of SDSS J1256— 0224 are consistent with 
membership in the Galactic inner halo, with estimated UVW space velocities indicating a slightly 
prograde, eccentric and inclined Galactic orbit (3.5 ^ ^ 11 kpc; \Zmax\ — 7.5 kpc). Comparison to 
synthetic spectra computed with the Phoenix code, including the recent implementation of kinetic 
condensate formation (Drift-Phoenix), indicate Te// « 2100-2500 K and [M/H] w -1.5 to -1.0 
for log 5 « 5.0-5.5 (cgs), although there arc clear discrepancies between model and observed spectra 
particularly in the red optical region. The stronger metal-oxide bands present in the Drift-Phoenix 
model spectra, a result of phase-non-equilibrium abundances of grain species, appears to contradict 
prior suggestions that grain formation is inhibited in metal-poor atmospheres; conclusive statements 
on the metallicity dependence of grain formation efficiency are as yet premature. In addition, an 
apparent shift in the temperature scale of L subdwarfs relative to L dwarfs may obviate the need for 
modified grain chemistry to explain some of the former's unique spectral characteristics. 

Subject headings: stars: chemically peculiar — stars: individual ( SDSS J125637. 13—022452.4) — 
stars: low mass, brown dwarfs — subdwarfs 
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1. INTRODUCTION 

L subdwarfs are the lowest-luminosity and least- 
massive halo popula tion dwarf s t ars c urrently known 
(jBurgasser. Kirkpatrick. fc Lipinel l2005l ). They derive 
their name from their gross spectral similarities to the lo- 
cal L dwarf population of very low mass stars and brown 
dwarfs fsee iKirkpatrickl 120051 ) . but are distinguished by 
specific spectral anomalies, including the presence of en- 
hanced metal hydride absorption bands and unusually 
blue near-infrared colors. These features are indicative of 
subsolar atmospheric abundances, as similar peculiarities 
distinguish metal-poor M subdwarfs from M dwarfs (e.g ., 
iMould fc H^da^[l97i [Gi^fT997^ et al.ll200oh . 

L subdwarfs, like M subdwarfs, exhibit kinematics consis- 
tent with membership in the Galactic halo, with inclined 
and eccentric Galactic orbits indicating dynamical heat- 
ing over long timescales and/or formation outside the 
Galactic disk (iDahn etld] 120081: IBurgasser eFall 120081 : 
iCushing et al.l [20091 ). The low luminosities (logLftoi/L© 

< —3) and effective tempe ratures of L subdwarfs (Te f f 

< 3000 K; iLeggett et atl I2OO0I: iReiners fc Basdl 12^61: 



IBurgasser et al.l [20081 ) . coupled with their non-solar at- 
mospheric abundances, are of particular interest for stud- 
ies of low temperature atmospheres, testing thermochem- 
istry and condensate format ion processes in chemically 
peculiar environments (e.g., lAckerman fc MarlevI I2OOII 
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lLodderd[200l iHelhng fc Woitkd[200l . In addition, as 
their inferred masses extend down to and below the 
metallicity-dependent hydrogen burning mini mum mass 
([Burrows et al.l 119931 : IBurgasser eFal] l2003a[ ). L subd- 
warfs can potentially test metallicity dependencies on 
low-mass star formation and brown dwarf evolution, rel- 
evant for instance in tracing the terminus of the m ain 
sequence in globular clusters (e.g.. [Richer et al.[|2008f) . 

Despite their utility to atmospheric, star for- 
mation and Galactic population studies, only 
four L subdwarfs have been reported to date: 

the prototyp e ^2MASS J05325346-H8246465 

([Burgasser et al.l l2003ai here after 2MASS 05 32-^8246), 
2MASS J16262034-h3925190 ([Burgassej|200l . hereafter 
2MASS J1626+3925 ), 2MASS J06164006-6407194 
(ICushing et all [2009L hereafter 2MASS J0616-6407), 
and SDSS J125637.13-022452.4 (iSivarani et all 
[2001 hereafter SDSS J1256-0224). The first 
three sources were identified s erendipitously i n the 
Two Micron All Sky Survey ([Skrutskie et al.l [20061 . 
hereafter 2MASS), and have been studied exten- 
sively at optica l and near-infrared wavelengths (e.g. 
Burgasser et al.l [20081: [Burgasser. Cruz fc Kirkpatrickl 



20071: ICush ing fc Vacca^ '2006': 'Gushi ng et all l200a 



Gizis fc Har vin 2006; Rcincrs fc Basri 20 06l:lPatten et al 
2006t IScholz et all 120091 : iSchilbach. Roser fc Scholz 
20M ). SDSS J1256-0224 was found in the Sloan 



Digital Sky Survey ([York et all [2000l hereafter SDSS) 
as pa rt of a directe d searc h for unusual red so urces . 
Both ISivarani et all ([2009f ) and IScholz et al] ([20091 ) 

have identified SDSS J1256-0224 as an L subdwarf 
based on its optical spectrum (L-type, with unusually 
strong metal hydride bands), blue near- infrared colors 
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{J ~ K, = 0.10±0.03; iSchilbach. R5ser fc Schoia [20091) . 
and halo-like kinematics. As yet, no detailed study of 
the optical and near-infrared spectral properties of this 
unusual source have been made. 

In this article, we present new observations of 
SDSS J 1256— 0224 and conduct a detailed analysis 
of its spectral and kinematic properties. Spectro- 
scopic observations spanning the red optical and near- 
infrared are described in § 2. The empirical prop- 
erties of SDSS J1256— 0224 are assessed in § 3, in- 
cluding classification (on the preli minary scheme of 
iBurgasser. Cruz fc Kirkpatricld l2007f ) , distance estima- 
tion, kinematics and Galactic orbit. In § 4 we ex- 
amine the atmospheric properties of this source by 
comparing its colors and spectra to the latest gener- 
ation of the Cond-Phoenix atm osphere simulations 
(jHauschildt. Baron fc AllardI [19971) and to the Drift- 
Phoenix model atm ospheres ()Dehn|[2007[: iHelling et all 
I2008bt [Witti 120081), the latter of which includes a ki- 
netic approach of phase-non-equilibrium dust formation. 
Results are summarized in § 5. 

2. OBSERVATIONS 
2.1. Red Optical Spectroscopy 

Optical spectra of SDSS J1256— 0224 were obtained 
on 2006 May 7 (UT) using the Low Dispersion Survey 
Spectrograph (LDSS-3) mounted on the Magellan 6.5m 
Clay Telescope. LDSS-3 is an ima ging spectrograph, up- 
grade d from the original LDSS-2 (jAUington-Smith et al.l 
|1994[ ) for improved red sensitivity. Conditions during 
the observations were clear with excellent seeing (0'.'6 at 
i'-band). The VP H- red grism (660 lines/mm) with a 
0'.'75 wide (4 pixels) longslit mask was used, with the slit 
aligned to the parallactic angle. This configuration pro- 
vides 6050-10500 A spectra across the entire chip with 
an average resolution of A/AA « 1800 and dispersion 
along the chip of ~1.2 A/pixel. The OG590 longpass 
filter was used to eliminate second order light shortward 
of 6000 A. Two slow-read exposures of 750 s each were 
obtained at an airmass of 1.12. We also observed the 
G2 V star G 104-335 {V = 11.7) immediately after the 
SDSS J1256— 0224 observation and at a similar airmass 
(1.14) for telluric absorption co rrection. The flux stan- 
dard LTT 7987 (a.k.a. GJ 2147: lHamuv"eraII[l99l was 
observed on the same night using an identical slit and 
grism combination. All spectral observations were ac- 
companied by HeNeAr arc lamp and flat-field quartz 
lamp exposures for dispersion and pixel response cali- 
bration. 

LDSS-3 d ata were reduced in the IRAF^ environment 
(fTbdv[[l986l) . Raw images were first corrected for am- 
plifier bias voltage, stitched together, and subtracted by 
a median-combined set of slow-read bias frames taken 
during the afternoon. These processed images were then 
divided by a median-combined, bias-subtracted and nor- 
malized set of flat field frames. The LTT 7987 and G 
104-335 spectra were optimally extracted first using the 
APALL task with background subtraction. The spec- 
trum of SDSS J1256— 0224 was then extracted using the 

* IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



G star dispersion trace as a template. Dispersion solu- 
tions were determined from arc lamp spectra extracted 
using the same dispersion trace; solutions were accurate 
to ~0.08 pixels, or ~0.1 A. Flux calibration (instru- 
mental response correction) was determined using the 
tasks STANDARD and SENSFUNC with observations 
of LTT 7987, which we have found provide sufficient cal- 
ibration to <10% over the 6000-9000 A spectral band 
([Burgasser. Cruz fc Kirkpatrickl [2007[ ). Corrections to 
telluric O2 (6855-6955 A B-band, 7580-7740 A A-band) 
and H2O (7160-7340 A, 8125-8350 A, 9270-9680 A) ab- 
sorption bands were determined by linearly interpolating 
over these features in the G dwarf spectrum, dividing 
by the uncorrected spectrum, and multiplying the result 
with the spectrum of SDSS J1256-0224. The two spec- 
tra of SDSS J1256— 0224 were then coadded to improve 
signal-to- noise, which ranged from ^15 at the 6600 A 
peak to a maximum of -45 at 8500A. 

The reduced red optical spectrum of SDSS J1256-0224 
is shown in Figure [U compare d to equivalent dat a for 
the sdM9.5 SSSPM J101 3-1356 (iScholz et al.l[2004a[) and 
2MA SS J1626-H3925 ijBurgasser. Cruz fc KirkpatTica 
[2001 .^ Our data for SDSS J1256-0224 have consid- 
erably higher sig nal-to-noise than the original SDSS dis- 
covery spectrum ([Sivarani et al.|[2009l ) and hi gher resolu- 
tion than contemporaneous observations by Schol z et al.l 
(2009). As originally pointed out by Sivaran fet al.l 
(|2009l) . the optical spectrum of SDSS J1256-0224 ex- 
hibits several characteristics indicative of an L dwarf, in- 
cluding an overall red spectral slope from 6000-8500 A; 
strong molecular bands of CrH (8600 A) and FeH (8700 
and 9900 A); and alkah fine absorption from K I (7700 A 
doublet), Na I (8182/8193 A doublet), Rb I (7798 and 
7946 A) and Cs I (8519 and possibly 8941 A). Equiva- 
lent width (EW) measurements for these lines are listed 
in Table [2] The K I doublet is extremely broadened, pro- 
ducing a V-shape notch in the spectrum that spans 7300- 
8100 A, also characteristic of L dwarf spectra. There 
are a number of peculiar features in the spectrum of 
SDSS J1256— 0224 that are not common to L dwarf spec- 
tra, however, including unusually strong bands of CaH 
(6900 A) and TiO (7200 and 8400 A), and numerous 
metal lines from Ca I (6571 A), Ca II (8541 A) and Ti I 
(7204, 8433 and 9600-9700 A^). The absence of these 
species in L dwarf spectra is largely attributed to the 
formation of Ca-T i and Ca-Al mineral con densates (e.g^ 
Allard et al.l[Mnl lLodde7^[200l [HeiliMTWoitke fc Thil 



20081) ■ and their presence in L subdwarf spectra has 



been interprete d as an indication of inhibited conden- 
sate formation ([Burgasser et al.ir2003at [Reiners fc Basril 
[2006( 1. Whether or not this is an accurate interpreta- 
tion (see § 4.3), unsually strong CaH and metal- line ab- 
sorption is a chara cteristic trait of inetal-poor M sub- 
dwarf spectra (e.g., [Mould fc Hvland [19761 : [Gizisl[1997f l 
and the presence of these features in the spectrum of 

^ These data were obtain ed with the Gemini Multi-Object Spec- 
trometer Pook et al.ll2004I V 

^ Absorption in the 9600-9700 A spectrum of 
2MASS J 1626-t-3925 was incorrectly associate d with TiH by 
IBurgasserl (2004f) . [Gushing & Vacca (2006) and IReiners fc Basril 
II2OO6I) have since identified these features as arising from the a 
^F-z ^F° multiplet of Ti I, and we adopt these identifications 
here. 
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SDSS J1256— 0224 supports its characterization as a 
metal-poor, low-temperature dwarf. 

Specific comparison of SDSS J1256-0224 to 
2MASS J1626-f3925 reveals remarkable similarities 
between these two sources, although the latter exhibits 
somewhat stronger FeH, CrH, Rb I and Cs I absorption 
features and somewhat weaker Na I lines. Indeed, 
the variation in features between the three spectra 
shown in Figure [1] suggests a sequence of very late-type, 
metal-poor subdwarfs, with SDSS J1256— 0224 having 
intermediate line and band strengths. The specific 
spectral classification of SDSS J1256— 0224 is discussed 
in detail in § 3.1. Note that no significant Ha emission 
or absorption is detected in any of these spectra. 

2.2. Near- Infrared Spectroscopy 

Low resolution near-infrared spectral data for 
SDSS J1256— 0224 were obtained in clear conditions 
on 2005 March 23 (UT) using the SpeX spectrograph 
(|Ravner et al.l l2003') mounted on the 3m NASA Infrared 
Telescope Facility (IRTF). We used the prism-dispersed 
mode of SpeX with a 0'.'5 slit (aligned to the parallactic 
angle), providing 0.75-2.5 /im spectroscopy with reso- 
lution A/ A A ~ 120 and dispersion across the chip of 
20-30 A pixel-i. SDSS J1256-0224 was observed at 
an airmass of 1.08. Four exposures of 180 s each were 
obtained in an ABBA dither pattern along the slit. 
The AO V star HD 111744 was observed immediately 
before SDSS J1256-0224 at a similar airmass (1.07) for 
telluric absorption and flux calibration. Internal flat 
field and Ar arc lamps were observed with the target 
and calibrator source for pixel response and wavelength 
calibration. 

Data w ere reduced using the SpeXtool package, ver- 
sion 3.1 (jCushing. Vacca. fc Rav'^ l2004( l using stan- 
dard settings. Raw science images were first corrected 
for linearity, pair-wise subtracted, and divided by the 
corresponding median-combined flat field image. Spec- 
tra were optimally extracted using the default set- 
tings for aperture and background source regions, and 
wavelength calibration was determined from arc lamp 
and sky emission lines. The multiple spectral obser- 
vations were then median-combined after scaling indi- 
vidual spectra to match the highest signal-to-noise ob- 
servation. Telluric and instrumental response correc- 
tions for the scien ce data were deter mined using the 
method outhned in IVacca et al.l ([2003), with line shape 
kernels derived from the arc lines. Adjustments were 
made to the telluric spectra to compensate for differ- 
ing H I line strengths and pseudo-velocity shifts. Fi- 
nal calibration was made by multiplying the spectrum 
of SDSS J1256-0224 by the telluric correction spec- 
trum, typically accu rate to within 10% across the 0.8- 
2.5 M m window (see iBurgasser. Kirkpatrick fc Burrow"i 
I2006f) . Instrumental response was determined through 
the ratio of the observed AO V spectrum to a scaled, 
shifted and deconvolved Kurucz^ model spectrum of 
Vega. Signal-to-noise ranged from ^80 at the J-band 
peak (~1 /im) to ~15 at i^-band (~2.2 /im). 

The reduced near-infrared spectrum of 
SDSS J1256— 0224 is shown in Figure^ again compared 
to equivalent SpeX prism data for SSSPM J1013-1356 

^ See |http : //kurucz ■ harvard ■ edu/ stars ■ html [ 



and 2MASS J1626-h3925 (|BurgasseH [200l . All three 
spectra show similar overall spectral morphologies, with 
strong molecular absorption features and relatively blue 
1.0-2.5 ^m spectral slopes. Again, 2MASS J1626+3925 
appears to be most similar to SDSS J1256— 0224 in 
terms of detailed spectral features. The strong FeH 
band at 0.99 /zm present in the optical spectrum of 
SDSS J1256— 0224 is clearly discerned in these data, 
intermediate in strength between SSSPM J1013-1356 
and 2MASS J1626-I-3925 but considerably stronger than 
any normal L dwarf. Strong FeH absorption is also 
present in the 1.2-1.3 fim region. H2O absorption is 
quite pronounced at 1.4 /im, but the 1.8 /um band is 
notably weaker. CO absorption at 2.3 /xm, a hallmark 
of M and L dwarf near-infrared spectra, is absent in 
the three subdwarf spectra shown here. The weakened 
1.8 /im H2O and 2.3 /im CO bands, and the overall 
blue spectral slope, can be attributed to enhanced H2 
absorption in all three sources, peaking at 2.1 / zm but 
spanning much of the 1.5-2.5 /xm regi on shown (Linsk^ 
1969 : Sa umon et al.lll994l : iBorvsowTJOr g ensen, fc Zheng 
|1997() . This absorption produces the very blue near- 
infrared color of SDSS J1256-022 4, J - = 0.10±0.03 
(jSchilbach. Roser fc Scholj l2009f )^. quite distinct from 
the colors of normal L dwarfs (J — Kg « 1.5-2.5, 
iKirkpatrick et al"]l2000( ). Atomic lines arising from the 
neutral alkalis Na I and K I are also seen, but the 
resolution of the SpeX prism data is insufficient to 
obtain meaningful EW measurements. 

3. CHARACTERIZATION OF SDSS J1256-0224 

3.1. Spectral Classification 

The most basic characterization of a late-type dwarf 
like SDSS J1256— 0224 is its spectral classification. 
However, while a well-defined red o ptical classifica- 
tion scheme exist for L dwa rfs (Kirkpa trick et aH 
Il999bl 120001 see also iGeballe et al. 2002 for discus- 
sion on the near-infrared classification of L dwarfs), 
there are simply too few L subdwarfs to define 
a robust scheme. We therefore followed the ap- 

Siroach outlined in Burgasscr. Cruz fc Kirkpatrickl 
20071 ), comparing the 7300-9000 A spectrum of 
SDSS J1256— 0224 to equivalent-resolution spec- 
tra of the L dwarf s pectral standards defined in 
IKirkpatrick et all (|1999bf ). Figure [3] shows the two best- 
matching standards, the L3 2MASS Jl 1463449-1-2230527 
and the L4 2MASS J l 1550087-^2307058 (data from 
IKirkpatrick et al.lll999bP 'l. Focusing on the 7300-9000 A 
region, the spectrum of SDSS J1256— 0224 appears to 
lie intermediate between these two standards, based on 
the height of the 7300 A spectral peak; the depth and 
breadth of the 7700 A K I doublet; and the depts of the 
8500 A TiO, 8600 A CrH and 8700 A FeH bands. There 
are discrepancies between SDSS J1256— 0224 and the 
L dwarf standards in this region, notably the spectral 
slope between 7800 A and 8400 A, the depth of the 
8200 A Na I doublet (stronger in SDSS J1256-0224), 
and the presence of additional Ca I, Ca II and Ti I lines 

^ Note tha t the 2MASS J - Ks = 0.66 color reported by 
ISivarani et"all l|2009h is in fact an upper limit; this source was 
not detected by 2MASS in the Ks band. 

^ These dat a were obtained with the Low Resolution Imaging 
Spectrograph l lOke et al.lll995l) . 
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in the spectrum of SDSS J1256— 0224. However, these 
differences are not nearly as extreme as those at shorter 
(i.e., the CaH and TiO bands between 6600-7300 A) 
or longer wavelengths (i.e., the strong H2 absorp- 
tion). The comparisons shown in Figure [3] indicate a 
spectral type of sdL3.5 for SPSS J1 256-0224 on the 
iBurgasser. Cruz fc Kirkpatrickl (|2007D scheme. This is 
only 0.5 subtypes earlier than the sdL4 classification of 
2MASS J1626+3925, consistent with both the overaU 
similarities between the spectra of these two sources 
and slight differences in their atomic line and molecular 
band strengths. 

3.2. Absolute Brightness and Distance 

Parallax distance measurements have recently become 
available for low-tem perature subdw arfs sp anning 
types sdM7 to sdL7 jMonet et al.l [1992; Dah n et al 



2004 IBurgasser eFan 120081 : ISchilbach. Roser fc Scholz 
20091) ■ including SPSS J125 6-0224 for which 
Schilbach. Roser fc Scholj (|2009D determine d — 



90±23 pc. This measurement is fairly uncertain, so 
we have compared it to the linear absolute magni- 
tude/spectral type relations recently quantified by 
[Cushing et al.. (2009) for ultracool subdwarfs: 



Mj = 8.02 + 0.313 X SpT 
Mh = 7.77 + 0.300 X SpT 
Mk, = 7.44 + 0.320 x SpT 



(1) 
(2) 
(3) 

where SpT(sdM7) = 7, SpT(sdLO) = 10, etc. These re- 
lation predict Mj = 12.23±0.22, Mr = 11.81±0.23 and 
Mk, = 11.75±0.25 for a type of sdL3.5±0.5, and hence d 
= 66±9 pc for SPSS J1256-0224 b ased o n the photom- 
etry of ISchilbach. R5ser fc Schold ()2009l ) (the distance 
uncertainty includes contributions from photometry, a 
0.5 subtype classification uncertainty, covariance matrix 
elements for the Cushmg et al. 2009 relations, and vari- 
ation in distance estimates between JHKg values) . This 
distance inferred esti mate is considerably close r than 
the mean value from ISchilbach. Roser fc SchoIS ()2009[ ) 
but nevertheless consistent within experimental uncer- 
ta inties. It is also roughly half the 120 pc estimate 
of lBurgasser. Cruz fc Kirkpatrickl (l2007l) an d 50% larger 
tha n the 42 pc estimate of ISivarani et al.l ()2009l ). As 
the ICushing et al.l (j2009r ) relations provide the most ac- 
curate assessment of the ultracool subdwarf spectral 
type/absolute magnitude scale thus far, we use our esti- 
mated value from these relations for subsequent analysis. 

3.3. Kinematics 

The p roper motion of SPSS J1 256-0224 as mea- 
sured bv ISchilbach. Roser fc Scholz] (|2009f) . m conjunc- 
tion with our estimated distance, yields a tangential 
velocity Vtan = 186±26 km s~^, where the uncer- 
tainty is dominated by the distance estimate. Note 
that this value is again roughl y half that estimated by 
iBur gasser. Cruz fc KirkpatricS (|2007 ,) due to the reduc- 
tion in the estimated distance. Nevertheless, this motion 
is still indicative of halo kinematics. A radial velocity 
(K) for SPSS J1256-0224 was computed from its opti- 
cal spectrum, by comparing the measured line centers of 
the atomic lines listed in Table |2| to vacuum wavelengths 
obtained from the National Institute of Standards and 



Tech nologies atomic line database^° ()Ralchenko et al.l 
I2008D . A hehocentric Poppler shift of -130±11 km s'^ 
was determined, which includes a bary centric motion 
correction of —15 km s~^. The uncertainty in Vr in- 
cludes the standard deviation in the line centers and a 
systematic uncertainty of 4 km s~^ based on the un- 
certainty of the wavelength dispersion solution. This 
value is nominally consistent wit h the — 90ib40 km s~^ 
cross-correlation measurement of ISivarani et al.l ([2009); 
our improvement in precision is likely due to our higher 
signal-to-noise spectral data. 

Adding in our radial velocity measurements, we 
find space velocities [U,V,W] = [-115±11, -101±18, 
— 150±9] km s""'^ in the Local Standard of Rest (LSR), 
assuming an LSR solar motion of \U, V, W](^ = [10, 5.25, 
7.17] km s~i (Pchnen fc Binnevi ll998l ). Note that we 
adopt a right-handed velocity coordinate system with 
positive U pointing radially inward toward the Galac- 
tic center. The large space motions in all three LSR 
components are again strong indications of Galactic halo 
membership for SPSS J1256-0224. 

3.4. Galactic Orbit 

To explore the kinematics of SPSS J 1256— 0224 is more 
detail, we calculated its Galactic orbit using the inferred 
UVW velocities as initial conditions. These velocities 
were first converted to a Gal actic inertial frame (assurn - 
ing Vlsr = +220 km s'^; iKerr fc Lvnden-Belll [l986h . 
and the position and distance of SPSS J1256— 0224 
transformed to galactocentric rectangular coordinates 
[X, Y, Z] aligned with [U, V, W], assuming a Solar po- 
sition of [-8.5, 0.027] kpc (IKerr fc Lvnden-Bel]lll986l : 
IChen et al.ll200ll ). We adopt the convention of positive 
X pointing toward the Galactic center to align with our 
definition for U. The Galaxy was modeled using a set of 
static potentials comprising a spherically-symmetric halo 
and bulge and an axisymmetric, thin exponential disk. 
This model is a simplified version of the one described in 
^ehncn fc Binncy (1998) and includes the three density 
distributions: 



Pbnlgc{r) ^ Pbo ( — 



Phalo{r) = PhO — 

\ah 

Pdisk(i?,2) = Soe-^/«^(5(z) 



1 



ah 



with spherical and polar coordinates r = 
VX^ + + and R = VX^ + Y^. With these 
simplifications, the potentials corresponding to the 
bulge and halo densities can be expressed in terms of 
special functions, and the potential of the disk evaluated 
on a grid by numerical integration. The parameters 
Pbo, ab, ab, Pho, ah, ah, Ph, So, Rd we re adopted as given 
for Models I and II in Table 2.3 of iBinnev fc Tremaind 
(2008), which fit the measured rotation curve of the 
Galaxy but bracket the range of allowable disk/halo 
mass ratio in the Solar circle (Model I describes a 
Galaxy dominated by the disk mass at the solar circle, 
in Model II the halo mass dominates at the solar circle) . 
The orbit of SPSS J1256-0224 was integrated using a 

|http: //physics .nlst . gov/asdS] 
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second-order leapfrog method (kick-drift-kick) with a 
constant timestep of 1 kyr for a total simulation time 
of 1 Gyr centered on the present epoch. Energy was 
conserved to better than 1 part in 10""* over the full 
length of the simulation, with the error dominated 
by the resolution of the grid on which the disk force 
and potential were interpolated. The Z component of 
angular momentum was conserved to 1 part in 10""'^'^. 

Figure d] displays the resulting orbit of 
SDSS J1256-0224 using Model I (Model II pro- 
vides essentially identical results). The general 
character of this orbit is similar to several other ultra- 
cool subdwarfs, with a prograde eccentric orbit with 
apoaps near the Solar radius (3.5 ^ R ^ 11 kpc, e 
= 0.5) and substantial deviations from the Galactic 
plane (Z^ax ~ ±7.5 kpc). In terms of eccentricity, 
the orbit of SDSS J1256-0224 is more similar to 
that of 2MASS J053 2-h8246 (e = 0.5, 3 < R < 8.5; 
iBurgasser et al.l 120081 ) and less ballistic than that of 
the sdM8 LSR 1425-^7102 (jLepine. Shara. fc RichI 
l2003f ) which plunges to within 1 kpc from the Galac- 
tic center (jPahn et al.l l2008l ). The inclination of 
SDSS J1256-0224's orbit, tanz « Rmax/Zmax, is quite 
a bit larger, nearly 45° with respect to the Galactic plane 
as compared to ~15° and ~30° for 2MASS J0532+8246 
and LSR 1425-1-7102, respectively. 

The orbital characteristics of SDSS J1256-0224, 
2MASS J0532-h8246 and LSR 1425-^7102 are aU con- 
sistent with membership in the Galaxy's imier hal o pop- 
ulation (jChibal^eeri [Mil; [GmHIq!^^ The 
inner halo is believed to dominate the halo population 
in the inner 10-15 kpc of the Galaxy and has a t ypi- 

cal metallicity o f [M/H] 1.6 (jCarollo et al.ll2007l : see 

also iGizii 119971 ). Membership in the inner halo sug- 
gests an orig in in the dissipative mergers of satellite 
galaxies (e.g., iSearle fc Zinn 1978; Chiba fc Beers 2003; 
IBell et all I2008D . It is possible that the similarity of 
these orbits arise from selection effects. Stars spend a 
larger percentage of their orbital periods near apoaps (in 
this case near the Sun), and short, highly eccentric or- 
bits would more frequently align with the Sun's Galactic 
position. However, the recent discovery of the L subd- 
warf 2MASS J0616-6406, whose highly retrograde orbit 
extends out to >30 kpc making i t a likely mernber o f 
the Galaxy's outer halo population (jCushing et al.|[2"009( ). 
suggests that ultracool subdwarfs may nevertheless be 
well-mixed in the vicinity of the Sun. 



4.1. Spectral Models 

Phoenix is a general-purpose model atmosphere 
code, using plane-parallel geometry, thermochemical 
equilibrium calculations and opacity sampling to self- 
consistently solve for the temperature-pressure profile, 
chemical abundances and radiative/convective energy 
transfer through the atmosphere. Several implementa- 
tions of Phoenix have been used to study late-type 
dwarfs, incorporating various assumptions on elemen- 
tal abundances, atomic and molecular opacities, line 
profiles, condensate formati on and c onvec ti ye overshoot 
(e.g. Hauschildt, AUard fc BaronI Il999t lAllard et all 
120011: Dohnas et al.ll2008l : iHelhng et al.ll2008b^ 

Here we examine two iniplement ations of Pho enix, 
the GAIA-COND m odels ( iHauschil dt et~al] l2003l) and 
the Drift models ()Dehnl 12007 : .Helhng et al.l l2008bt 
IWittd |2008[ ). which differ only in their treatment of 
the dust cloud layers. The former are most sim- 
ilar to the Cond-Phoenix model set developed by 
lAllard et al.l (|200lD . in which condensate species are 
treated as element sinks only and phase-equilibrium is 
assumed. Hence, Cond-Phoenix models simulate dust- 
free but element-depleted atmospheres. The Drift- 
Phoenix models apply an advanced model of non- 
equilibrium grain formation, including seed formation, 
growth, evaporation, sedimentation and convective up- 
mixing to simulate the size distribution, abundances and 
vertical distribution of grains and their material compo- 
sition IWoitke fc Helling 2003, 2004; Helling fc Woitki 
120061: iHelling. Woitke fc Thil |2008| ). In this approach, 
each size-variable grain is made of a variety of com- 
pounds which changes with height according to its 
formation history. Alternate prescriptions for mod- 
eling condensate grain formation in cool dwarf at- 
mosph e res have bee n expl o red bv lAckerman fc MarlevI 
((2f)0l'): lAl iard et al.l (l2003l): iCoop er et al.' C2003D: Pfsuiil 
[2002, 200 51): iTsuii. Nakaiima fcY anagisawa [^W^ and 



[Burrows. Sudarskv fc Hubenvl ( 20061) : a t horough com- 
parison of these cloud models is given in IHelling et al.l 
(2008a). 

For this study, we employed both GAIA Cond- 
Phoenix and Drift-Phoenix atmosphere models span- 
ning 2000 < Te// < 3500 K (steps of 100 K), logg = 5.0 
and 5.5 (cgs) and -3.0 < [M/H] < 0.0 (steps of 0.5 dex; 
Solar [M/H] = 0). E lemental ab undances are scaled from 
lAnders fc Grevessi (jl989i) and iGrevesse. Noels fc Sauvall 
( 1991) ! 



4. ATMOSPHERIC PROPERTIES 

The distance and kinematics of SDSS J1256-0224 
do not provide useful constraints on its atmospheric 
properties — T^eff, surface gravity (logg) and metallic- 
ity ([M/H]). Such determinations require empirical cali- 
brations (e.g., well -characterized coeva l companions or 
cluster properties; iGizis fc ReidI Il997f) or direct com- 
paris on to spectral models (e.g.. iWoolf fc W allersteinI 
120061 ). As SDSS J1256-0224 is a seemingly isolated 
source, we used the latter approach, employing the 
most recent generat ion of the Cond-Phoe nix atmo- 
sphere simulations (IHauschildt. Baron fc AU ard 1997; 
IHauschildt. Allard fc BaronI Il999t iBaraffe et al.i i2003i) 
and the Drift-Phoenix model atm ospheres ([Dehnl 
[20071 : IHelling et al.l[20f)8bl : [Witt3[200l . 



4.2. Color Comparisons 

We pursued a two-step comparison of 
SDSS J1256-0224 to the models, first examin- 
ing optical/ near- infrared colors i' — J and J — 
(see also IScholz et al.l l2004bt iDahn et all 120081 : 
iSchilbach. Roser fc Schold 120091 ). Synthetic colors 
were computed directly from the model atmosphere 
spectra by convolving each with filter profiles from SDSS 
([Fukugita et al.l 119961 on the AB magnitude system) 
and 2MASS ([Cohen. Wheaton fc MegeathI 120031 on the 
Vega magnitude system). The filter profiles include 
the effects of detector quantum efficiency, telescope 
throughput and atmospheric transmission, all of which 
are e ssential for the complex sp ectra of late-type dwarfs 



are e ssential lor tne complex sp ' 
(;e.g.. lStephens fc Legge"ttil2004[ ) 
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Figure [5] displays colors for both model sets for log g 
= 5.5 and the full range of Te// and [M/H]. In addition 
to measurements for SDSS J1256— 0224, we all show 
color data for five subdwarfs cla ssified sdM8 and later: 
the sdM8 LSR 1425+7102, (jLepine. Shara. fc RichI 
[20?)^: the sdM8.5 2MAS S J 01423 153+0523285 
(|Burgasser. Cruz fc Kirkpatrickl l2007l here- 
after 2MASS J0142+0523); SSSPM J1013-1356, 
2MASS J1626+3925 and 2MASS J0532+8246. Near- 
infrared photometry for al l sourc es are from 2MASS or 
iSchilbach. R5ser fc Schoi3 (120091 ) . with the exception 
of 2MASS J0142+0523 for which a synthetic J - 
color was calculated using near-infrared spectral data 
from iBureasser et all (l200l . SDSS i' magnitudes for 
SDSS J1256-0224 and 2MASS J1626+3925 are from 
SDSS DR6; for the other sources, i' magnitudes were 
bootstrapped from In phot ometry from t he Su per- 
Cosmos Sky Survey (SSS; iHamblv et all [2001al rblB 
by calculating synthetic i' — In colors from pubhshed 
optica l spectral data for the se sources (Burgasser ct alJ 
2003at iL^pine. Shara. fc Ric h 2003; Scholz ct al. 2004b,; 



Burgasser. Cruz fc Kirkpatri ck 2007). The uncertainties 



for the spectroscopy-based magnitudes were assumed to 
be 0. 1 mag, based on prior work (e.g., iBurgasser et al.l 
[200^ . 

For the GAIA Cond-Phoenix models, predicted col- 
ors encompass the measured values of all the subdwarfs 
shown with the exception of 2MASS J0532-H8246. Lower 
temperatures for a given metallicity generally yield red- 
der i' — J and bluer J ~ Kg colors, except at for [M/H] 
< -2.0 for which i' — J colors actually turn blue. Lower 
metallicities at a given Te/ / lead to bluer i' —J and J— Kg 
colors. SDSS J1256-0224 and 2MASS J1626-I-3925 both 
fall along the [M/H] = -1.0 line in this diagram (as do 
LSR 14254-7102 and 2MASS J0142-h0523) around Teff 
~ 2300 and 2150 K, respectively. For log 5 — 5.0, colors 
for these sources agree with model atmosphere metal- 
licities and temperatures ~0.2 dex and '--^lOO K lower, 
respectively. The inferred model parameters based on 
these colors must be treated with caution, however, as 
the [M/H] = GAIA Cond-Phoenix models do not 



track well with mean 



J versus J — colors of MO- 



LD field dwarfs as compiled bv IWest et all (|2008l ). The 
divergence is likely due to absence of condensate opacity 
in the GAIA models, which are known to be necessary in 
reproducing the co lors of late- type M and L dwarfs (e.g., 
lAllard et al.l[200lh . 

The Drift-Phoenix models exhibit very different 
color trends below Tg/ / w 2500 K as condensates become 
a prominent opacity source in the photosphere, both 
i' — J and J — Kg colors trend redder for lower Te/ / for 
[M/H] > -2.5, with a notable kink in color tracks at Te// 
« 2200-2300 not present in the GAIA Cond-Phoenix 
models. The additional reddening places the [M/H] = 
models into closer agreement with the mean colors 
of LO dwarfs, although they still diverge from M5-M9 
dwarfs. The J— Kg color reversal effects all metallicities, 
with the result that none of the models reproduce the 
measured colors of the L subdwarfs SDSS J1256— 0224, 
2MASS J1626-h3925 and 2MASS J0532-)-8246 (log 5 = 
5.0 models show similar behavior). The deviation of the 
models away from the measured photometry is likely due 
to the model i-band magnitudes, which are highly sen- 
sitive to the strong molecular opacity present at these 



wavelengths (see § 4.3). A rough extrapolation of the 
models suggests lower metallicities for SDSS J1256— 0224 
and 2MASS J1626-I-3925, [M/H] - -2.0 to -1.5. 

4.3. Spectral Comparisons 

To further assess the agreement between models with 
observations, we compared the observed spectrum of 
SDSS J1256-0224 directly to the GAIA Cond-Phoenix 
and Drift-Phoenix model spectra for the same range 
of parameters shown in Figure O These comparisons 
were made to the combined red optical and near-infrared 
spectrum, which was stitched together by first smoothing 
the individual spectra to a common resolution of A/ A A 
= 100, scaling the spectra to match flux densities in the 
0.8-0.9 /im range, and then combining the red optical 
spectral data for A < 0.9 /im with the near-infrared data 
for A > 0.9 fim^^ (see Figure [2|). The uncertainty spec- 
trum (flux uncertainty as a function of wavelength) was 
combined using the same scaling factors and wavelength 
cutoffs. We then interpolated the entire spectrum onto 
a linear wavelength scale (note that a wavelength scale 
uniform in frequency produced similar results). 

Observational and model spectra were initially normal- 
ized over the 0.9-1.0 /im range; then, for each model 
spectrum, a goodness-of-fit statistic was calculated, 

[F(A)-a%|(A)]2 



^^'^ ^ a%}(A)a(A) 



(4) 



Here, F(X) and (t(A) are the observed spectrum and un- 
certainty; S'{p}(A) is the model spectrum for model pa- 
rameters {p} = {Te//, logg, [M/H]}; a is a normaliza- 
tion scale factor for the model spectrum; and W{X) is a 
weighting function that satisfies J2{\}WW — 1- This 
form was chosen as a compromise between a standard 
formulation (e.g., J2[-^ ~ ^]'^/^) ^^'^ ^ reduced for- 
mulation (e.g., ~ 'S']^/"'^)? as the former places too 
much emphasis on the lowest signal regions (i.e., strong 
absorption features) while the latter places too much em- 
phasis on the highest signal-to-noise continuum regions 
(e.g., the 0.9- 1.1 At m peak) . For al t ernat e approaches, see 
iTakedal (Il995h and lCushing et all (l2008h . The sums were 
computed over the spectral range {A} = 0.64-2.4 fj,m. 
The normalization factor a was allowed to vary over 0.5- 
1.5 to account for continuum offsets between the observed 
and model spectra, and the normalization with the min- 
imum r was retained. Various weighting functions H^(A) 
were considered, but ultimately we settled on one that 
was constant for all wavelengths.^^ Note that we do not 
consider P a robust estimator; it merely provides a quan- 
titative measure of the best-fit model to the data. The 
best fits (minimum F) were also visually compared to 
verify that they did indeed provide a good match to the 
data. 

Table [3] lists the parameters for the five best fits for 
the GAIA Cond-Phoenix and Drift-Phoenix mod- 
els; the single best-fit models are compared with the data 

Red optical data were only used up to 0.9 /^m due to concerns 
over the relative flux ca libration of these data over the 0.9—1.0 fiia 
range; sec discussion in Bur gasser. Cruz &: KirkpatricH 11200 71) . 

Gushing et al. (2008), in their analysis of optical and infrared 
spectra of late-type dwarfs, considered a weighting function that 
scaled with the width of each spectral wavelength bin. As we in- 
terpolate the observed and model spectra onto a common, linear 
wavelength scale, our constant weighting scheme is equivalent. 
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in Figure [6] The best-fit model parameters are similar 
between the two sets, with Teff — 2300-2500 K and 
[M/H] = -1.5 to -1.0 for the GAIA Cond-Phoenix mod- 
els, and Teff = 2100-2400 K and [M/H] = -1.5 and for 
the Drift-Phoenix models. These parameters includes 
models with logg = 5.0 and 5.5; lower gravities are gen- 
erally matched with lower metallicities, which may sim- 
ply indicate a tradeoff in the inferred photospheric pres- 
sure (Pph oc. g/K, where k, is the Rosseland mean opac- 
ity which is generally smaller for lower metallicities). A 
higher surface gravity is in fact preferred if this object 
is a low-mass member of the halo population as evolu- 
tionary models predict that a 5 Gyr source with Te/ / = 
2100-2500 K should have log . <? = 5.3-5.5 and mass 0.08- 
0.085 Mq (jBurrows et al.ll2001l : iBaraffe et al.ll2003D . In- 
deed, there are more best-fitting models with log 5 = 5.5 
than 5.0. 

Encouragingly, the best-fit parameters for both GAIA 
Cond-Phoenix and Drift-Phoenix models are sim- 
ilar to the best-fit parameters from the GAIA Cond- 
Phoenix color comparison. The metallicities inferred 
from both color and spectral comparisons are also con- 
sistent with the me an metallicities of inner halo stars 
(|Carollo et al.|[2007[) . In addition, we find that both sets 
of models do a reasonably good job at matching the over- 
all spectral energy distribution of SDSS J1256— 0224, in 
particular fitting the blue spectral slope from 1.3-2.4 /im 
and the depth of the 1.4 and 1.9 /im H2O bands. The 
GAIA Cond-Phoenix models also provide a fairly good 
match to the 0.9-1.3 /im spectral peak, reproducing the 
strong FeH bands at 0.99 and 1.25 /im but predicting 
excessively strong alkali lines in this region. The Drift- 
Phoenix models do a poorer job in this region, failing to 
reproduce the 1.1 /xm spectral peak and, like the GAIA 
Cond-Phoenix models, exhibiting excessively strong al- 
kali lines. 

The are more significant deficiencies in the red opti- 
cal region, however, with both model sets failing to re- 
produce spectral features in detail, particularly around 
the spectral peaks at 6600 and 7400 A. The models 
predict excessively strong Rb I and Na I alkali lines, 
and appear to be missing CrH and FeH opacity in 
the 8600-8700 A region. The Drift-Phoenix mod- 
els exhibit excessively strong TiO absorption at 8400 A, 
and pronounced discrepancies in the 6700 A CaH and 
7200 A TiO bands. Surprisingly, the older GAIA Cond- 
Phoenix models provide overall better fits to the red 
optical data, although there are still clear problems in 
alkali line strengths, excessive emission in the pseudo- 
continua between 6000-7500 A, and missing CrH and 
FeH opacity. These discrepancies, which influence model 
i' — J color trends (Figure O are probably attributable 
in part to inadequate treatment of the far- wing line pro- 
files of the pressure-broadened Na I and K I doublets, 
as more recent opacity cal culations have not incorpo- 
rated in the present models ([Burrows fc Volobuvevl[200l 
iJohnas et al.l l2008f ) . The treatment of dust formation 
can influence these line proflles due to fee dback on the 
temperature structure ([Johnas et al.ll2008l ). It is there- 
fore promising that the K I line cores and other al- 
kali line strengths are reproduced better in the Drift- 
Phoenix models. Incompleteness in molecular opaci- 
ties (e.g., TiO, CrH and FeH) are also likely responsi- 
ble, a well-known problem in modeling the optical spec- 



tra of normal L dwarfs (e.g., iKirkpatrick et aLlll999af) . 
We cannot rule out the additional influence of elemen- 
tal composition variat ions, which are present in subsolar 
meta l licity stars (e.g . . lEdvardsson et al.lll993l : iFulbrigh^ 
[20001 : lAsDlund et al.! Vmi^ and can influence the over- 
all atmospheric chemical pathways. The red optical 
region clearly remains problematic for low-temperature 
model atmospheres regardless of m etallicity (see also 
iBurgasser. Cruz fc Kirkpatricld l2007l ). although the in- 
corporation of condensate dust formation appears to be 
aiding alkali line fits somewhat. 

4.4. Discussion 

While it appears that the inclusion of condensate grain 
formation as specified by the Drift-Phoenix models 
provides some improvement to the near-infrared col- 
ors and alkali line profiles of ultracool dwarfs and sub- 
dwarfs, the overall better spectral and color fits to 
SDSS J1256-0224 by the condensate-free GAIA Cond- 
Phoenix models suggests that grain chemistry may be 
unimportant in metal-poor atmospheres. This has been 
the conclusion of several studies, citing the presence of 
enhanced metal-oxide bands, strong lines from refractory 
species, and blue near-infrared colors as consistent with 
largely condensate - free photo s pheres ([Burgasser et al.l 
2003a; .Gizis fc Ha7^ l2006t iReiners fc Basril l2006t 
IBurgasser. Cruz fc Kirkpatricldl2007[ ). 

However, there are problems with this simple in- 
terpretation. It is clear that TiO and VO features 
weaken from the M subdwarfs to the L subdwarfs, 
and continue to weake n even through the current L 
subdwarf sequence (e.g.. IBurgasser. Cruz fc KirkpatricS 
l2007f ). This trend is consistent with the depletion of 
refractory gas-phase elements with decreasing temper- 
ature. In contrast, TiO and VO bands are stronger 
in later-type M giants wh ere low atmosp heric pressure 
inhibits condensation (see lLoddersll200"4 l. Conclusions 
that gaseous TiO and VO bands are nevertheless en- 
hanced may also be biases by si mple equilibrium treat- 
ments of condensation chemistry. iHellin g. Woitke fc Thil 
(2008) have shown that rare-element compounds, includ- 
ing Ca- and Ti-bearing condensates, never achieve phase- 
equilibrium in low-temperature atmospheres; hence, the 
abundances of gas molecules are generally higher in non- 
equihbrium cloud models such as Drift-Phoenix. The 
fact that the metal-oxide bands observed in the spec- 
trum of SDSS J1256— 0224 are actually weaker than 
predicted by the best-fit Drift-Phoenix cloud mod- 
els suggests that condensate grain formation may actu- 
ally be more efficient in metal-poor atmospheres than 
these models predict. There are some important caveats 
to this interpretation, however. Incomplete molecular 
and pressure-broadened line opacities are clearly an issue 
for low-temperature atmospheres, affecting in particular 
pressure-temperature profiles and associated chemistry 
through the photosphere. Also, we have only examined 
one prescription of grain formation in this study; cur- 
rent cloud models h ave not yet reached agreement (e.g., 
iHelling et al.l [2008a[ ). It is clear that conclusive state- 
ments on the efficiency of condensation in metal-poor 
atmospheres are premature; models are simply not yet 
adequate to address this question. 

Another consideration, independent of the state of 
current model atmospheres, is the possibility that L 
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subdwarfs are simply warmer than equivalently-classifed 
L dwarfs. The Te// inferred for SDSS J1256-0224, 
roughly in the range 2100-2500 K, is co mparable to 
those for solar- metallicity L0-L2 dwarfs (e.g. JVrba et al.l 
|2004[ ). At these temperatures, condensates play 
a less prominent role in atmospheric opacity than 
for cooler/later- type dwarfs. Similarly, the sdL7 
2MASS J0532+8246 has been shown to have a Te// = 
1730±90 K, comparable to L4-L5 dwarfs. This shift is 
due in la rge part to the clas sification methodology pro- 
posed by iBurgasser. Cruz fc Kirkpatr ick (2007). Com- 
parison of L dwarf and subdwarf optical spectra in the 
7300-9000 A range emphasizes the importance of the 
pressure-broadened 7700 A K I doublet. This feature 
is inherently pressure-sensitive, and at a given tempera- 
ture will be deeper and broader in the higher-pressure, 
metal-poor (small k) photospheres of L subdwarfs. As 
such, the apparent persistence of TiO absorption in 
L subdwarfs may simply reflect a shift in tempera- 
ture scale and not c ondensation abnormalities (see also 
IBurgasser fc Kirkpatr ick 2006). 

Such shifts in inferred quantities like temperature and 
metallicity have been seen as an inherent weakness in 
the classi fication scheme o f ultracool subdwarfs in gen- 
eral fe.g.. iJao et al.l (|2008f )). However, the current state 
of flux in theoretical models should also emphasize the 
importance of divorcing classification (a purely empiri- 
cal exercise) from interpretive physical parameters. The 
former is built upon a set of specific standard stars, not 
models, although standards are as yet in insufficient sup- 
ply for the L subdwarfs. As more examples are un- 
covered, measurement of their luminosities and Te//S , 
and improved theoretical modeling particularly of opti- 
cal spectra, will provide necessary the constraints to cali- 
brate a future classification sequence and enable more ro- 
bust assessment on condensate grain formation efficiency 
that is currently possible. 

5. SUMMARY 

We have presented a thorough analysis of the red 
optical and near-infrared spectrum of t he L subdwarf 
SPSS J1256-0224, originally identified bv lSivarani et all 
(|2009l) in the SDSS survey. This source is simi- 
lar to the sdL4 2MASS J1626-f3925 at both optical 
and near-infrared wavelengths, and we determine an 
sdL3.5 classification following the pre liminary scheme of 
IBurgasser. Cruz fc Kirkpatrick! ()2007l ). Using the abso- 
lute magnitude/spectral t ype relations for ultra cool sub- 
dwarfs recently defined bv lCushing et al.l ()2009[ ). we esti- 
mate a distance of 66±9 pc to this source, formally con- 
sistent with the less p recise 90±23 pc parallax distanc e 
measurement made bv lSchilbach. R5ser fc Schold (|2009D . 



Combined with its high proper motion and radial veloc- 
ity, we confirm that SDSS J1256— 0224 is a kinematic 
member of the Galactic inner halo, with a modestly ec- 
centric and highly inclined Galactic orbit whose apoaps 
is near the Sun. A comparison of the colors and spec- 
tra of SDSS J1256-0224 to GAIA Cond-Phoenix and 
Drift-Phoenix atmospheric models indicate best-fit at- 
mospheric parameters of T^// = 2100-2500 K and [M/H] 
= -1.5 to -1.0 for log5 = 5.0-5.5, although discrepancies 
between the models and the data, particularly in the red 
optical region, indicate that these parameters be treated 
with caution. Comparisons to the Drift-Phoenix mod- 
els contradict prior conclusions that condensate forma- 
tion may be inhibited in metal-poor, low-temperature 
atmospheres. Non-equilibrium grain species abundances 
(particularly for Ti- and Ca-bearing species) predict even 
stronger metal-oxide bands than those observed in the 
spectrum of SDSS J1256— 0224, and thus enhanced con- 
densate formation in the atmosphere of this metal-poor 
source. Improvements to the model atmospheres in the 
red optical region are necessary before any conclusive 
statement can be made. In addition, the possibility that 
the temperature scale of L subdwarfs is warmer than 
that of L dwarfs may provide a sufficient explanation for 
strong metal-oxide bands irrespective of grain chemistry. 
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Fig. 1.— Red optical spectra of the SSSPM J1013-1356 (top), SDSS J1256-0224 (middle) and 2MASS J1626+3925 (bottom). Spectra 
are normalized in the 8500-8600 A region and offset for comparison (dotted lines). Primary spectral features are indicated, as well as 
regions of strong telluric absorption (e) in the spectra of SSSPM J1013-1356 and 2MASS J1626+3925. 
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Fig. 2.— Low resolution near-infrared spectra of SSSPM J1013-1356 (top), SDSS J1256-0224 (middle) and 2MASS J1626+3925 (bottom), 
all obtained with SpeX in prism mode. Joined to these are the optical data from Figurc[T] smoothed to the same resolution (A/AA ~ 120). 
Spectra are normalized in the 0.9-1.0 fim region and offset for comparison (dotted lines). Primary spectral features are indicated. 



L Subdwarf SDSS J125637.13-022452.4 



13 




6500 7000 7500 8000 8500 9000 

Wavelength (Angstroms) 

Fig. 3.— Optical spectrum of SDSS J1256-0224 (black lines) compared to the L dwarf optical standards (red lines) 2MASS J1146+2230 
(L3) and 2MASS J1155+2307 (L4). All spectra arc normalized in the 8250-8350 A and offset for clarity (dotted hnes). In the 7300-9000 A 
range, the spectral morpholog y of SDSS J1 256— 0224 is intermediate between the two standards, indicating an sdL3.5 spectral type for this 
source based on the scheme of IBurgasser. Cruz fc Kirkpatrick (,2007,). 
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Fig. 4. — Galactic orbit of SDSS J1256— 0224 over 1 Gyr centered on the current epoch, based on Galactic Model I from lBinnev fc Tremairi^ 
II2008I ). Upper left and right pan els shows orbit i n \X, Y] and \R, Z ] inertial frame coordinates; the current position of the Sun (Xq = 
-8.5 kpc, Zq = +27 pc; IKerr fc L vnden-Bell 198|; |Ch en et al.1l2001l) is indicated by the symbol. Bottom panels show time evolution of 
R and Z. In all panels, past motion is indicated by dashed lines, future motion by solid lines, and current position by the black point. 
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Fig. 5.— SDSS-2MASS i' - J versus J - Ks colors for the late-type subdwarfs 2MASS J0142+0523, LSR 1424+7102, SSSPM J1013- 
1356, SDSS J1256-0224, 2MASS J1626+3925 and 2MASS J0532+8246 as compared to GAIA Cond-Phoenix (top) and Drift-Phoenix 
(bottom) atmospheric model predictions. Models are shown for logg = 5.5, 2000 < T^eff 5: 3500 K in steps of 100 K (along solid lines) 
and -3 .0 < [M/H] < in steps of 0.5 dex (along dotted lines). Also shown are mean SDSS-2MASS colors of MO-LO dwarfs from lWest eTaTI 
ll2008lV 
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Fig. 6. — (Left) Best-fit GAIA Cond-Phoenix (top) and Drift-Phoenix models (bottom; red lines) compared with the combined 
red optical and near-infrared spectrum of SDSS J1256— 0224 (black lines). Both models and spectra have been smoothed to a resolution 
A/AA = 100 and interpolated onto a common wavelength scale. Observational data are normalized in 0.9—1.0 /im range while models are 
normalized so as to minimize residuals. Best-fit parameters are listed in Table |3] (Right) Comparison of best-fit models (based on full 
spectrum) to observed data in the red optical. Here, both models and spectra have been smoothed to a resolution A/AA = 1000 and 
interpolated onto a common wavelength scale. The goodness-of-fit statistic F computed over the 6400—9000 A range is notably worse than 
that for the full optical/ncar-infrared spectrum. 
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TABLE 1 
Properties of 
SDSS J125637. 13-022452. 4. 



PtircLrri6ter 




Ref. 


"J2000 


12'»56™3?;16 


1 


l5j2000 


-02°24'52'.'2 


1 


Spectral Type 
i' -J 


sdL3.5 


2 


3.253±0.024 


3,4 


J-Ks 


0.097±0.026 


4 




12.23±0.22 


2,5 


Mh" 


11.81±0.23 


2,5 




11.75±0.25 


2,5 


dest (pc)'' 


66±9 


2,4,5 


Vtan (km S~^) 


186±26 


2,4 


Vr (km s"i) 


-130±11 


2 


U (km s-l) 


-115±11 


2 


V (km s-i) 


-101±18 


2 


(km s-i) 


-150±9 


2 


Te// (K) 


~ 2100-2500 


2 


logg (cgs) 


- 5.0-5.5 


2 


[M/H] (dex) 


~ -1.5 - -1.0 


2 



Refere nces. — (1) 

2MASS (iSk rutskie et~all [200g i: 

2) This paper; (3) SPSS 

Adelmari - McCarthy ot al. 2008|); 
4) ISchilbach. Ro scr & Sc holj 

20091 ): (5) iCushing eTa l. (2009). 
^ Estimated absolute magnitudes 
based on the absolute magni- 
tude/spectral t ype relations of 
ICushing et aD l|2OO90 and spec- 
tral type sdLS.S.'' Note that 
[Schilbach, Roscr & Scholz (2001) 
measure an astrometric distance of 
90±23 pc for this source, formally 
consistent with our ^more pr ecise 

J stima te based on the'Cushing ct aLl 
20091 ') absolute magnitude/spectral 
type relations. 



TABLE 2 
Atomic Line Equivalent Widths 
(EW). 



Feature 


Lino Center 


EW 




(A) 


(A) 


Ha 




>-0.9 


Ca I 


6571.10 


3.1±0.6 


Til 


7204.47 


2.4±0.5 


Rb I 


7798.41 


2.9±0.4 


Rb I 


7945.83 


2.6±0.2 


Na I 


8182.03 


9.1±0.2'=- 


Na I 


8193.33 




Til 


8433.31 


3.7±0.3 


Cs I 


8519.19 


1.3±0.2 


Ca II 


8540.39 


0.5±0.3 


Cs I 


8941.47 


1.1±0.4 



^ EW for combined doublet. 
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TABLE 3 
Spectral Model Fits. 



Model 


log 9 
(cgs) 


(K) 


[M/H] 
(dex) 


F 


GAIA Cond-Phoenix 


5.5 


2300 


-1.0 


0.51 




5.5 


2400 


-1.0 


0.53 




5.0 


2400 


-1.5 


0.56 




5.0 


2300 


-1.5 


0.57 




5.5 


2500 


-1.0 


0.72 


Drift-Phoenix 


5.0 


2300 


-1.5 


0.64 




5.0 


2200 


-1.5 


0.66 




5.5 


2300 


-1.5 


0.67 




5.5 


2400 


-1.5 


0.75 




5.5 


2100 


-1.5 


0.76 



Note. — Top five best fit models for 
combined optical and near-infrared spectrum of 
SDSS J1256-0224, ranked by F. Best-fit models 
are shown in Figure [§] 



